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Abstract 

The ocean thermohaline circulation, also called meridional overturning circula- 
tion, is caused by water density contrasts. This circulation has large capacity of 
carrying heat around the globe and it thus affects the energy budget and further 
affects the climate. We consider a thermohaline circulation model in the meridional 
plane under external wind forcing. We show that, when there is no wind forcing, the 
stream function and the density fluctuation (under appropriate metrics) tend to zero 
exponentially fast as time goes to infinity. With rapidly oscillating wind forcing, we 
obtain an averaging principle for the thermohaline circulation model. This averaging 
principle provides convergence results and comparison estimates between the origi- 
nal thermohaline circulation and the averaged thermohaline circulation, where the 
wind forcing is replaced by its time average. This establishes the validity for using 
the averaged thermohaline circulation model for numerical simulations at long time 
scales. 

Mathematics Subject Classifications: Primary 35K35, 60H15, 76U05; Sec- 
ondary 86A05, 34D35 

Key Words: Exponential decay, averaging principle, geophysical flows, wind forcing 

1 Introduction 

In addition to the wind-driven surface circulation, the ocean also exhibits a large merid- 
ional overturning circulation cahed the thermohaline circulation. The ocean is heated 
(thus made less dense) where pure freshwater is evaporated (water thus made saltier 
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and denser), and vice versa. The global thermohaline circulation involves water masses 
sinking at high latitudes and upwelling at lower latitudes. The process is maintained 
by water density contrasts in the ocean, which themselves are created by atmospheric 
forcing , namely, heat and water exchange via evaporation and condensation. During the 
thermohaline circulation, water masses carry heat (or cold) around the globe. Thus, it 
is believed that the global ocean thermohaline circulation plays an important role in the 



climate [28|. 



A two-dimensional thermohaline circulation model involves the Navier-Stokes equations 
for momentum (in the meridional plane) together with the convection-diffusion equations 
for temperature and for salinity. Due to the linear equation of state (relating fluid 
density with temperature and salinity), these latter two equations may be replaced by a 
single convection-diffusion equation for density or density fluctuation. We consider the 
thermohaline circulation on a fluid domain in the vertical meridional xz-plane, as in, e.g., 

A^/;* + J(AV', V) = Px + i^^'^ip + f{x,z,t) , (1.1) 
pt + J{p,^) = -N^^, + ^Ap, (1.2) 

where il^{x, z, t) is the stream function which defines the velocity field (u, w) = {ipz, — V'lOi 
p{x, z, t) is the density fluctuation from the mean density; > is the viscosity; A^^ > 
is the mean buoyancy frequency and is taken as a constant; Pr is the Prandtl number. 
Finally, J(a, b) = Cxby — aybx is the Jacobian operator and A = dxx + dzz is the Laplace 
operator. Note that ^ + J{-,'4') = ^ + '^^x + wdz is the material derivative. The wind 
forcing term f(x,z,t) is to be specified below. 

In [Q], the author made some numerical simulation for (1.1), (1-2) with periodic boundary 
conditions in both x and z for ip and p and / = 0. Therefore, we assume that ■i/' and P 
are periodic (with period 1) in x and z, and also assume that ip and p have zero mean. 



In some recent work on the thermohaline circulation, the wind forcing is ignored |27, 



29, In this paper, we will consider the impact of wind forcing on the thermohaline 
circulation, while considering the evolution of the fluid density fluctuation (rather than 
the fluid density itself). 

In the first part of this paper, we obtain the exponential decay estimates for the stream 
function ip and density fluctuation p with / = (no external wind forcing). In the 
second part, we consider the effect of wind forcing on the stream function and density 
fluctuation. We obtain an averaging principle for rapidly oscillating wind forcing, which 
provides convergence results and comparison estimates between the original thermohaline 
circulation and the averaged thermohaline circulation. This establishes the validity for 
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using the averaged thermohahne circulation model for numerical simulations at long time 
scales. 



2 Exponential Decay: Without Wind Forcing 

In this section, we consider the long time behavior of the stream function and the den- 
sity fluctuation in the thermohaline circulation. We first briefly comment on the local 
existence for (1.1)-(1.2) with periodic boundary conditions (with period 1) in both x and 
z for ij) and p (with zero mean), we introduce some notations: 

dxdz^ 

D 

where D = {(x, z) : < x, z < 1} is the periodic fluid domain. 

H = LU = {n: n € L^(i?),n is perodic both in x and z,/n = 0}, with norm || • ||; 

V = iJpg^o = : u£ H,\/u£ H}, H^^^Q = {u : u £ H,Vu £ H,Au £ H}, etc. 
In fact, by a result in Q, we know ||Au|| is equivalent to 

Define the vorticity u! = Aip. It is well known that A^""^ exists for A with periodic 
boundary and zero mean, then (1.1) can be written as 

LUt + J{u;, A~^a;) = Px + vAu + /(x, z, t). (2.1) 

Since the nonlinear Jacobian term is continuous from V xV ^ H x H^hy the theory of 
[11|, we have the following local existence result for (2.1) and (1.2): 

Lemma 2.1 (Xiocal Existence^ Let {ujq, pq) £ V x V (initial values, that is tpo £ H^^^q) 
and f £ L°°{0,T; H), then (2.1) and (1.2) with periodic boundary conditions in both x 
and z for ip and p with zero mean has a unique local solution satisfying 

u £ L~(0, T; V) n L2(o, T; HI^^^) p £ L°°(0, T; V) H L°°(0, T; H^^,^), 
that is, (1.1) and (1.2) has a unique local solution satisfying 

£ L°°(0, T; Hl,,^) n ^^(o, T; F^Vo) P e ^~(0, T; V) n L°°(0, T; F^^.q), 
where T depends on {ujq,pq). 

We need the following properties and estimates (see Q) of the Jacobian operator J : 
Hq X Hq ^ in the sequel: 

Jif,g)hdxdy = - J{f,h)gdxdy, J{f, g)gdxdy = 0, 
D Jd Jd 
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J{f,g)dxdy\ < IIV/IIIIV5II, 

D 

for all f,g,h£ Hq. 

Now, we derive some a priori estimates for the solution of (1.1) and (1.2) with / = to 
ensure the global existence. For / 7^ 0, the estimates are almost the same. 

Multiplying (1.1) by ip, performing an integration by parts and using the periodic bound- 
ary conditions, we conclude 

^|||VV|P+ |5.pV + ^^I|AV|P = 0. (2.2) 

Similarly, we get 

Multiplying ( |2^ ) by iV^ and adding to (|2^), we have 

1 |(iV2|| + \\pf) + v{N^\\Mf + -^llVpf ) = 0. (2.4) 

By the Poincare inequality and taking a = ^min{l, -p^KAi is the smallest eigenvalue 
of —A with periodic boundary and zero mean), we obtain 

^'llVVf + \\P\? + j\{N^\\^i:f + ^^\Wpf)]dt 



2 dt '""'I ~^ J ^""P"^ + -p-r^\^P\\ = 0- (2.3) 



<e-*(iV^||VVor + ||por). (2.5) 

If we only want to know whether the solution tends to zero as t — > 00, we could use the 
following special Gronwall Lemma. We omit the details for this asymptotics here, and 
we will concentrate on the exponential decay of the solution in the sequel. 

Lemma 2.2 // a non-negative differential function f satisfies 

f'{t) + aif{t)<g{t), 
where ai > and lim g{t) = 0, then lim f{t) = 0. 



Proof. We first have 



f{t) < e~"i*/(0) + e-"i* r e'''^g{T)dT. (2.6) 

Jo 
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Since lim g{t) = 0, for any given e > 0, there exists a T such that 

t~»oo 

g{t) < e, for t > T. 

So, (2.6) can be written as 

f{t) < e-"i*/(0) + — (1 - e-"i*) + e-"i* e"i^g(T)dT. (2.7) 

«! JO 

Let t — > oo and note that e is arbitrary, the result is thus obtained. 
Now, we turn to the exponential decay estimates for the solutions. 

Multiplying (1.1) and ( |1.2| ) by —Ap and Atp respectively, integrating by parts, using the 



periodic boundary conditions, we have 

1 d 



1 d 

Note that 



--IIAV'f = J d^pA^p-i.\\VAi;f, (2. 
~\\Vpf-Jjip,i;)Ap = -N'J d,i;Ap-^\\Apf. (2. 



Atpf = -J pAiP^ - uWVAi, 



ld_ 

2di 

1 



< - (1 -5)zy||VAV'f ( VO < 5 < 1) 

1 „ „2 (1 -(J)z^„ . ,„2 

where the Poincare inequality is used. By the Gronwall inequality and (2.5), we thus 
obtain 

<1-S)iy, X, nt (1-5)1/ 



~ U{l-5)u^ Jo -""^(^^llWor + IIPoine^'dr. (2.10) 

Noticing that a = min{l, -p^}, we know for every given Pr, there exists a. < 5 < 1, 
such that 

Hence (2.10) can be reduced to 
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II AV^f < e-^^dlAV'of + 4^(r3|j^(^'l|VV'of + llPof )). (2.11) 

Multiplying (2.8) by N'^, adding to (2.9), using integrating by parts and the fact that 
/[/ dxV pV'ip + / dxVil^V p\ = 0, we conclude that 

i|(iV2||AV'f + llVpf ) - / J(p, V)Ap = -iVVllVAV'f - ;^||Apf . (2.12) 
Now we only need to estimate — / J(p, V') A/9 = / {pzipx — Pxi^z){Pxx + Pzz)- Note that 



J PzlpxPxx = ^ J pI^I-^xz - J PxPzi^xx, 
J Pzi^xPzz = J Pllpxz, 



Pxi^zPxx = ^ J Pli'xz, 
-J Pxi^zPzz = ~\j pI'^xz + J PxPzi^zz- 

Now wc need the following lemma about the equivalence between ||?/||ff2 and ||Ati|| for 
Lemma 2.3 For every u G -H^erO? ^^^^ 



||ti||j^2 < ai||Att||, 

where ai = ^1 + Ai + Af . 
Proof. Since 

Adding || filial to both side of (2.13), and using Poincare inequality, the proof of this 
lemma is complete. 

Using this lemma, we imply that 

I j J{p, ^)Ap| = \j{pl- pl)i>xz - j PxPzii'xx - Az)\ 

< I {\px? + \pz\'^Mxz\ + I IPo^II/O^KIV'xxI + \tpzz\) 

< 



a,{2 + V2)\\Aml |Px|')^+ / IP.I')^)- 
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By the following inequality from |^] 

1 1 1 \/2 1 

||m||l4 < a2||u||5||Vu||5,M G Hl^^Q{D),a2 = + — + 2)i 

Therefore, 

' J(ft V')Ap| < o;<!^(2 + V2)||Av||||Vp||||Ap|| 



Combining (2.12) with (2.14) and using the Poincare inequality, we obtain 

~{N'\\Mf + \Wpf) 
< -NM^Mf + ( "i«2(6 + 2V2)Pr _ s,)-^)\\Vpf. (2.15) 

By (2.11), there exists a ti such that 



(6 + 2^/2)Pr 



AV-f - (1 - S2)—fr = -<53 < 0, < ^3 < 1. 



4^21^ " " ' XiPr 

Let $0 = IIAV'of + 4^(1-1)3,.^ (^'ll^V'of + ll/^of )• It follows that 

_(!_£). 4(^2(1-52)1^^ 

e ^1 $0 < ^ 



a\a^{<o + 2^/2)\lPr'^' 
Since 4^2(1 — 82) < l, so ti can be chosen as 

Ai ^ .aafaf(6 + 2^/2)AiPr2$0i 
(1 — o)v 

By (2.15), we know 

||Vpf < Mo, < t < ti (2.16) 

and 

llVpf < e-"^*(iV2||^^^||2 ^ ||^^^||2)^ ^ > (2.17) 

Here 02 = minj^s, "^XTJ" ^'^'^ ^® constant depending on ti, ||A'i/'o|Pi z^,Pr and 

Ai. 

The estimates (2.16) and (2.17) tell us that the mean-square norm of the density (fluc- 
tuation) gradient, ||V/9||, is uniformly bounded up to some time instant and then decay 
exponentially fast. 
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Remark 2.4 // u is large enough or Pr is small enough, we also can have 

402 XiPr 

that is 



Pr < 



alal{2 + ^ li^t^n^^^H^ + + \\Po\?))) 

or V satisfying 

e-^^dlAV'of + i^(r^|)^(^'l|V^of + llPof ))) < 1 

i.e., 



V > 



laiaj{6 + 2^/2)Pr'^Xl 
4(52(1 -(52) ' 



Now, we derive the estimates for ||VAV'|| when the initial value V'o € -ffperO- -^y ^ similar 
process as above, we have 



llVAV-f - J J{M, = -J d^pA^iP - lyWA^iPf. (2.18) 



l_d 

By the same discussion as in deriving (2.14), we have 



ly J(t;,'ii)A?;| < a?a^(2 + v^)||A'u||||Vi;||||Ai;||, for u,v e H^^rO- (2-19) 
We use inequality (2.19) for: u = ^ and v = Aip, we have 

1 1 J(AV',V)A2VI<a?ai(2 + V2)||AV'||||VAV^||||A2v.|| 

< ^WA^^f + MMi^i^llAV'fllVAVf (2.20) 



and 



I \^,pA'i^\<^\\A'^pf + l\\Vpf. (2.21) 
By (2.18), (2.20) and (2.21), we get 



Thus we conclude by using the Poincare inequahty, 

ilVA^f < ^llVpf + (^M&tM) II A^f - ^)||VA^f . (2.22) 
at u V 2 

Using (2.11) and for Ai large enough, there exists as > such that 

I AV'II^ - All/ = -as < 0, (a2 > "3). 



aiai{12 + AV2)..^ ,,,2 



Hence 



|||VA^f + a3||VAV^f <^||Vpf, 



by the Gronwall's inequality, we have 

llVAV-f <e-"3*||vAV;of + - T e'^^^*"^) ||Vpf dr. 

Jo 

For t < ti, using (2.16), we have 



IIVAV-f < e-^'^tiiVAV'of + —Mo = Mi. (2.23) 

lyas 

For t > ti, using (2.17) and a2 > as, we get 

IIVAV'II^ < e-"3*||VAV'o||^ + -(iV^||AV'o||^ + ||Vpo||^) /* e-"3(*-^)-"2r^^ 

Jo 

= e-^'^^WVA^of + -{N^\\A^of + \\Vpof)) [ e-^"^^-"^^^^ dr. 

V Jo 

< e-"3*(||VAV^o||' + ^(iV2||AV^o||' + ||Vpo|P)). (2.24) 
If Ai is not so large, by (2.11), there exist 0:4 > and t2 > large enough such that 

afa|(12 + 4V2)„ , ,„o , ^ , 

^ ^ ^IjAV'lr - Aiz^ = -a4 < 0, for t > t2 (^2 > Q!4)- 

Thus (2.22) can be written as 

^llVAV^f + a4||VA^f <^||V/9f. 

By the the same argument as we obtain the estimates of (2.23) and (2.24), we could get 

similar estimates. 

The estimates (2.23) and (2.24) tell us that the mean-square norm of the vorticity gradi- 
ent, ||Vu;||, is uniformly bounded up to some time instant and then decay exponentially 
fast. 

Thus, by the above estimates (2.5), (2.11), (2.16), (2.17), (2.23) and (2.24), which hold 
in the case of no wind forcing (/ = 0), we obtain the main theorem in this section: 
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Theorem 2.5 (Exponential Decay in the Case of No Wind Forcing^ Let the 

initial conditions for the vorticity and density fluctuation (ujo,pq) be in V x V (i.e., 
ipo G r.). Then, when there is no external wind forcing, llipW iis and llplliifi tend 

^ perO perO 

to zero exponentially fast as t —> oo. That is, under appropriate norms or metrics, the 
stream function and density fluctuation tend to zero exponentially fast as time goes to 
infinity. 

Moreover, there exists some time instant T > such that, the mean-square norms for 
the density (fluctuation) gradient, ||Vp||, and for the vorticity gradient, \\\/uj\\, are uni- 
formly bounded when t < T and exponentially decay when t > T. The time instant 
T = max{ti,t2} depends on Pr, N"^, v, Ai and initial values. 

3 Averaging Principle: Rapidly Oscillating Forcing 

In this section, we consider the averaging principle for the system of (2.1) and (1.2) under 
the rapidly oscillating forcing f{x,z,t). We rewrite (2.1) and (1.2) as 

uJt-\-Aiuj = -J{uj,A^'^uj)-\-px-\-f{x,z,t), (3.1) 
Pt + A2P = -J{p,A~'^uj) - N^A'^uj^ , (3.2) 

where and A2 denote the operator — z^A and — -^A with the periodic boundary 
conditions and zero mean. For the rest of this section, we concentrate on the system 
(3.1)-(3.2). 

We assume that the forcing term / in (3.1) is rapidly oscillating, i.e., it has the form 
f{x, y, t) = /(x, y, rjt) = f{r}t), with parameter rj ^ 1. We also assume that / has a well- 
defined time average. With such a forcing, it is desirable to understand the fluid dynamics 
in some averaged sense, and compare the averaged flows with the original (non-averaged) 
flows. 

The main result of this section is an averaging principle for (3.1)-(3.2) with rapidly 
oscillating forcing on finite but large time intervals. This includes comparison estimate 
and convergence result (as rj — > 00) between (3.1)-(3.2) and its averaged motions. 
Starting from the fundamental work of Bogolyubov [Q the averaging theory for ODE 
has been developed and generalized in a large number of works (see 0-|9| and the 
references therein). Bogolyubov's main theorems have been generalized in |jl^ to the 
case of differential equations with bounded operator-valued coefficients. Some problems 
of averaging of differential equations with unbounded operator-valued coefficients have 
been considered in [11|-|14] in the framework of abstract parabolic equations. In |15], 



Ilyin considered the averaging principle for an equation of the form 
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dtU = N{u) + f{r]t), (3.3) 

where / is a given forcing function and ^ 1 is a large dimensionless parameter, and / 
has a time average defined as 

lim ^ f f{s)ds = fo. (3.4) 



t—*oo t Jo 

Note that and A2 are sectorial operators. For a sectorial operator, one can define the 



fractional power of A as follows |11]: 



1 f°° 

j^a ^ where ^"'^ = — - / e-^e'^^dt. 

r(a) Jo 

The corresponding domains D(A") are Banach spaces with norm given by 

\\x\\a ■■= \\x\\d{A'^) = 
We recall some definitions and results to be used in the rest of this section. 



Lemma 3.1 [H] The following estimates are valid: 

\\e--^^\\L2^L2 < Ke-^\ t > 0, (3.5) 

\\A^e-^'h2^L2<^e-'^\ t>0, (3.6) 
where K, Ka are positive constants. 



Remark 3.2 Since Ai and Ai are different operators but both satisfy the conditions of 
Lemma 3.1. For the simplicity, we take the same constants for Ai and A2 when we use 
Lemma 3.1. 



Lemma 3.3 |J1|/ Given two sectorial operators A and B in , let D{A) = D{B), 
Rea{A) > 0,Rea{B) > 0, and for some a G [0,1). Let the operator {A - B)A~°' 
be bounded in . Then for every 7 G [0,1), D{A^) = D{B'^), the two norms being 
equivalent. 

Setting 

T = r]t, e = r/"\ 
we rewrite the equations (3.1)-(3.2) in the so-called standard form 
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oJr + eAiUJ + eJ{A ^u>,ll)) = ep^ + ef{x,y,T), (3-7) 
Pr + eA2p + eJ{p,A-'^uj) = -eN'^A'^u^ , (3.8) 

We assume that / has a time average , fo{x, z), in D^A^); the value of 7 will be specified 
later on. More precisely, let /(t),/o G A"^ and suppose that 

\\A^^ 1'^^ fiT)dT - fo)\\ <mm{M^,a^{T)), (3.9) 

where > 0, cr^(T) ^0, as T — > 00. 
We consider the averaged equation 

u!r + eAiui + eJ{A^^u),u!) = epx + efo{x,y), (3.10) 
Pr + eA2p + eJ{p,A-^Ld) = -eiV^A^^o;^ , (3.11) 

By the method of |jl^-|l^], we know the semigroup St corresponding to equation (|3.10| )- 

1 1 

( |3.11| ) possesses absorbing sets in the space H = Lp^^Q x Lp^^Q,Y = D{Ai) x D{A'2) = 
H^^^Q X i?ig^o and D{A) = D{Ai) x D{A2){m fact, D{Ai) = D(^2))-Using Lemma 3.3, 
we know D{AJ) = D{A2), for 7 E [0, 1]. || • || and || • || 1 denote the norm in L^g^Q and 
iJpg^o- These sets are certain balls B{Ro) in these spaces, where Rq is large enough. This 
means that for every bounded set B 

StB C B{Rq), for t > to{B,Ro). 

In addition, the semigroup is uniformly bounded in these spaces, that is, given any ball, 
in particular, the ball B{Rq), there exists a ball B{R) such that 

StB{Ro) C B{R), for t > 0. 

By increasing R we may assume that 

StB{Ro) C B{R-r), fort > 0,r > 0, 

where r is a positive constant. We consider the averaging principle in the space V. 
Given a point luq in Bv{Rq), we compare the trajectories (solutions) {uj{t), p{t)) and 
(w(r),/)(r)) of system (3.1-(3.2) and (3.7)-(3.8) starting from same initial point. Con- 
sider their difference on the interval r G [0, > ^ being arbitrary but fixed. We suppose 
for the moment that {lj{t), p{t)) £ Bv{R). Then the difference z{t) = uj{t) — 0{t), Q = 
p — p satisfies the equations 

drZ + eAiz{T)+e[J{A-^u;,u;)- J{A-^iu,iu)] = e((p, -/>,) + (/(r) - /o))(3.12) 
dre + eA2{T)+e[J{A-'^iu,p)-J{A-^LU,p)] = - e { A''^ lo - A-'^lu, (3.13) 
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We first give the following lemma, the proof can be obtained by direct estimate. 

Lemma 3.4 The nonlinear operator J{u, v) is a bounded Lipschitz map in the following 
sense: 

\\J{ui,Vl) - Jiu2,V2)\\ < 

Cl(\\ui\\l + ||'U2|| 1 + Ikllli + Ik2|| l)(|kl — ■Willi + ||'"2 — V2\\l), (3-14) 
2 2 " " 2 " " 2 11112'^" "2 "2 

where Ci is some positive constants. 

2 

Inverting the linear operators Ai and A2 we come to the equivalent integral equations of 
(3.12) and (3.13) 

z{t) = -e r e-^-^^^^-'^[J{A-^uj,io)-J{A-'^Ld,Ld)]ds 
Jo 

+ e r e-'-^'^^-'\p, - -p.)ds + e r e~'-^'^^-'\fis) - fo)ds, (3.15) 
Jo Jo 

e(T) = -e r e-'-^'^^-'\j{A-^u,p)- J{A-^u,p)]ds 
Jo 

- eN^ r e-'^'^^-'\A-^u^- A-^u^). (3.16) 







Using (3.6) and (3.14), we see that the || • ||i-norm of the first term in the right hand 
side of (3.15) satisfies the inequality 

||e r Afe-^^^^^-'^[J{A-^io,Lo) - JiA-^ud,Cd)]ds\\ 
Jo 

< e r KiCie-^(T - syh-'^^^-'hrnizis^ids 
Jo 2 2 ' 2 

= 2RKiCie^ r (t - syh-^'^^-'hlzisMlids. (3.17) 

2 2 Jo 2 
Let us estimate the second term in the right hand side of (3.15). 

||e r Aie-'-^'^^-'\p^ - p^)ds\\ <Kie^ /^(r - s)-5e-^"(^-") ||e(s)|| 1 ds. (3.18) 

Jo 2 Jo 2 

Now let us estimate the third term in the right hand side of (3.15). Integrating by parts 
we have 

lie r e-^^^(--'\f{s)-fo)dsh_ 
Jo 2 

= II - ee-^-^^(--^) J\f{t) - fo)dt\l + £ Aie-'^'^^-'^ J\fis) - fo)ds\\ 1 

< ||e4"V^-^^M? r{f{t)-fo)dt\\ 
Jo 
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+ lle^ ^e-'^^^^-^Uj Hfis) - fo)ds\\. (3.19) 

Jo Js 

Using (3.6) and (3.9), we further have 

WeAr'^e-'^'^Aj Hfit) - fo)dt\\ < eKi e— -(er)^-5||i r AUfit) - fo)dt\\ 
Jo 2 ' r Jo 



< (er)3+'^Ki_^mm(M^,(j^(r))e-'"^ =: L(r). (3.20) 

For any 5 > 0, let ts be so large that for r > Ts,wehavea^ < 6. Let eo be so small that 
for e < Co the inequality ^ > is valid. Then 



L(r)<G^i(T,e) = e- 



(er)i+^Ki M^, iiT<Ts. 



Let 7 > — i. Note that ts does not depend on e. Taking S and then e — 0, we 



2 

obtain 



/"(/(*) - /o)dt||i < G^i(T,e) ^ when e ^ 0. (3.21) 

Jo 2 



Jo Js 2 ' Jo 

■2 ^ Jo 2 T Jo 

= Ki_^{-f + ^)-i(M^(er^)H7 + ^tH^) =: g^^{T, e), (3.22) 

where for any > wc have chosen ry;^ so large that cr-y(r) < fi when t > t^. Letting 
ji —> and then e ^ we obtain 

G^2(r,e)^0, e^O. 
Thus, by (3.17)-(3.22) we obtain the following inequality: 

||^(t)||i <i^e5 r {t - s)-"2{\\z{s)U + \\Q{s)h_)ds + G^{T,e), (3.23) 

2 Jo 2 2' 

where K = max{2i?KiCi } and = Gyi + G^2 ^ 0, e ^ 0. 



2 2 2 ' 
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Similar to the argument above for z{t), and using (3.6) and (3.14), we get 



|e||i < Kie^ Ht - s)-^i\\zis)\\i + \\eis)\\i_)ds, (3.24) 

2 JQ 2 2 



where Ki depends on R,Ki,Ci and Ai. 

2 2 

Adding (3.23) and (3.24), we finally obtain 







\\z{T)\\i+\\e{T)\\^<{K + K^)e-2 I {r-srH\\z{s)U + \\Q{s)U)ds + G,{T,e). (3.25) 
Here we need the following fact. 



Lemma 3.5 J7|/ Let 7 G (0, 1] and for t £ [0, T] 

u{t) <a + b [ {t- sy'^u{s)ds. 

Then 







u{t) < aE^i{bT{^))-t), 
where the function E^{z) is monotone increasing and E^{z) ~ j~^e^ as z ^ 00. 

Applying this lemma to the inequality (3.25) on r € [0, we obtain 
Mt)\\l + \mr)h- < 

2 2 

GJT,e)Ei{eTTT{K + Kif) < GJT,e)Ei{TTT{K + Kif):=r]l{e). (3.26) 

2 2 

We thus have proved the proximity of solutions of (3.1) and (3.2) in V, for the trajectory 
{uj{t),p{t)) with initial condition (u;(0),p(0)) G i?v(-Ro) stays in the ball B{R) on the 
time interval [0, ^1- 

Let e be so small that the right-hand side of (3.26) are less than |, where r is defined earlier 
in this section when we discuss absorbing sets. Suppose that the trajectory {u;{t),p{t)) 
leaves the ball B{R) during the interval [0, ^] and let r* be the first moment where 

ll'^(''"*)l|i + ll/'(''"*)ll 1 = However, on the interval r G [0,r*l both trajectories stay in 

2 2 

the ball B{R) and what we have proved so far shows that the inequality 

\\u:{r)-u:{r)h_ + \\p{r)--p{T)U<'- 
2 2 z 

is valid. In particular, it is valid for r = r*. This together with the inequality ||a;(T*)|| 1 + 

2 

IIp(t*)|| 1 < R—r, which holds by the hypothesis of the following theorem and the property 

2 

of the semigroup S{t), gives the contradiction 

Mr*)h + Mrnh- < ||^(r*)-c^(r*)|U+||p(T*)-/^(T*)||i+||cD(r*)||i+||/^(r*)|U </?-^, 
22 2 2222 
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since ||u;(r*)||i+||p(r*)||i =it;. 

Thus we have the main theorem in this section: 

Theorem 3.6 ('Averaging Principle in the Case of Rapidly Oscillating Wind 
Forcing^ Assume that the wind forcing has a time average. Let T > be arbitrary 
and fixed. If j > — ^ and the initial conditions for the vorticity and density fluctua- 
tion (uj{0), p{0)) = {lD{0), p{0)) are in the absorbing ball i?v(^o) (depending on then 
for T G [0,^], we have the following comparison and convergence estimate between the 
thermohaline circulation and the averaged thermohaline circulation 

Mt) - cD(t)||i + \\p{t) - p(r)|| 1 < ?7^(e) ^ 0, as e ^ 0, 

2 2 

where r}\,{e) is a decaying function defined in (3.26). 

4 Summary 

The ocean thermohahne circulation has important impacts on the cHmate. We have 
considered a thermohahne circulation model in the meridional plane under external wind 
forcing. We have shown that, when there is no wind forcing, the stream function and 
the density fluctuation (under appropriate metrics) tend to zero exponentially fast as 
time goes to infinity (Theorem 2.5). With rapidly oscillating wind forcing, we have 
obtained an averaging principle for the thermohaline circulation model (Theorem 3.6). 
This averaging principle provides convergence results and comparison estimates between 
the original thermohaline circulation and the averaged thermohaline circulation, where 
the wind forcing is replaced by its time average. 
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